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We consider the evolution of a flat Friedmann-Robertson- Walker Universe, filled with a causal 
bulk viscous cosmological fluid, in the presence of variable gravitational and cosmological constants. 
The basic equation for the Hubble parameter, generalizing the evolution equation in the case of con- 
' stant gravitational coupling and cosmological term is derived, under the supplementary assumption 

, that the total energy of the Universe is conserved. By assuming that the cosmological constant is 

' proportional to the square of the Hubble parameter and a power law dependence of the bulk vis- 

cosity coefficient, temperature and relaxation time on the energy density of the cosmological fluid, 
O , two classes of exact solutions of the field equations are obtained. In the first class of solutions the 

' Universe ends in an inflationary era, while in the second class of solutions the expansion of the 

Universe is noninflationary for all times. In both models the cosmological "constant" is a decreasing 
function of time, while the gravitational "constant" increases in the early period of evolution of the 
Universe, tending in the large time limit to a constant value. 
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I. INTRODUCTION 



Recent observations of type la supernovae with redshift up to about z < 1 provided evidence that we may hve in 
a low mass-density Universe, with the contribution of the non-relativistic matter (baryonic plus dark) to the total 
energy density of the Universe of the order of il„i ^ 0.3 0]- The value of il„i is significantly less than unity Q| 
and consequently either the Universe is open or there is some additional energy density p sufficient to reach the value 
^totai = 1, predicted by inflationary theory. Observations also show that the deceleration parameter of the Universe 
^1 ^ • g is in the range —1 < 9 < 0, and the present-day Universe undergoes an accelerated expansionary evolution. 

Several physical models have been proposed to give a consistent physical interpretation to these observational facts. 
• • ' One candidate, and maybe the most convincing one for the missing energy is vacuum energy density or cosmological 
.5^ ' constant A Q. 

. Since the pioneering work of Dirac j^, who proposed, motivated by the occurence of large numbers in Universe, a 
^ ' theory with a time variable gravitational coupling constant G, cosmological models with variable G and nonvanishing 
and variable cosmological term have been intensively investigated in the physical literature [Q- In the isotropic 



cosmological model of Chen and Wu |14| it is supposed, in the spirit of quantum cosmology, that the effective 
cosmological constant A varies as (with a the scale function). In the cosmological model of Lima and Maia |2^] 
the cosmological constant A = A [H) ~ ?>j3H^ + 3 (1 — /3) /Hi is a complicated function of the Hubble parameter 
a constant j3 and an arbitrary time scale HJ^ , leading to a cosmic history beginning from an instability of the de 
Sitter space-time. The cosmological implications of a time dependence of the cosmological of the form A ^ have 
been considered by Berman ||l2| . Waga investigated flat cosmological models with the cosmological term varying 
as A = a/ a? -f I3H^ -\- 7, with a, (3 and 7 constants. In this model exact expressions for observable quantities can 
be obtained. Nucleosynthesis in decaying-vacuum cosmological models based on the Chen-Wu ansatz has been 
investigated by Abdel- Rahman The consistency with the observed helium abundance and baryon asymmetry 
allows a maximum vacuum energy close to the radiation ener gy today. Anisotropic Bianchi type I cosmological 
models with variable G and A have been analyzed by Beesham and it was shown that in this case there are no 
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classical inflationary solutions of pure exponential form. Cosmological models with the gravitational and cosmological 
constants generalised as coupling scalars and with G ~ a" have been discussed by Sistero Generalized field 

equations with time dependent G and A have been proposed in and in an attempt to reconcile the large 
number hypothesis with Einstein's theory of gravitation. Limits on the variability of G using binary-pulsar data have 
been obtained by Damour, Gibbons and Taylor A detailed analysis of Friedmann- Robertson- Walker (FRW) 

Universes in a wide range of scalar-tensor theories of gravity has been performed by Barrow and Parsons ||l0|] . 

Dissipative effects, including both bulk and shear viscosity, are supposed to play a very important role in the early 
evolution of the Universe. The first attempts at creating a theory of relativistic fluids were those of Eckart |^ and 
Landau and Lifshitz |^ . These theories are now known to be pathological in several respects. Regardless of the choice 
of equation of state, all equilibrium states in these theories are unstable and in addition signals may be propagated 
through the fluid at velocities exceeding the speed of light. These problems arise due to the first order nature of the 
theory, that is, it considers only first-order deviations from the equilibrium leading to parabolic differential equations, 
hence to infinite speeds of propagation for heat fiow and viscosity, in contradiction with the principle of causality. 
Conventional theory is thus applicable only to phenomena which are quasi-stationary, i.e. slowly varying on space 
and time scales characterized by mean free path and mean collision time. 



A relativistic second-order theory was found by Israel |30| and developed by Israel and Stewart Hiscock and 



Lindblom ||32| and Hiscock and Salmonson |33| into what is called "transient" or "extended" irreversible thermody- 
namics. In this model deviations from equilibrium (bulk stress, heat flow and shear stress) are treated as independent 
dynamical variables, leading to a total of 14 dynamical fluid variables to be determined. For general reviews on causal 
thermodynamics and its role in relativity see |34j and . 

Causal bulk viscous thermodynamics has been extensively used for describing the dynamics and evolution of the 
early Universe or in an astrophysical context. But due to the complicated character of the evolution equations, very 
few exact cosmological solutions of the gravitational fleld equations are known in the framework of the full causal 
theory. For a homogeneous Universe filled with a full causal viscous fiuid source obeying the relation ^ ~ p^^^, with 
p the energy density of the cosmological fiuid, exact general solutions of the field equations have been obtained in 
p6| - [H - It has also been proposed that causal bulk viscous thermodynamics can model on a phenomenological level 
matter creation in the early Universe |39| , . Exact causal viscous cosmologies with ^ ~ p" ,s ^ 1/2 have been 
considered in |^ . 

Because of technical reasons, most investigations of dissipative causal cosmologies have assumed FRW symmetry 
(i.e. homogeneity and isotropy) or small perturbations around it |^5|| . The Einstein field equations for homogeneous 
models with dissipative fiuids can be decoupled and therefore are reduced to an autonomous system of first order 
ordinary differential equations, which can be analyzed qualitatively |^^, |^^, 48 1. 



Recent developments in particle physics and cosmology have shown that the cosmological constant A ought to be 
treated as a dynamical quantity rather than a simple constant. The dynamics of the scale factor in FRW type models 
with a variable A term has been recently revisited for the pefect fluid case by Overduin and Cooperstock |4^ . 

The effects of dissipation as expressed in the form of a non- vanishing bulk viscosity coefficient in the stress-en ergy 
tensor of the matter in cosmological models with variable A and G have been considered by several authors |Q- Ml - 
The role of a transient bulk viscosity in a FRW space-time with decaying vacuum has been discussed in ||5^ . Models 
with causal bulk viscous cosmological fiuid have been considered recently by Arbab and Beesham jH^ . They obtained 
both power-law and infiationary solutions, with the gravitational constant an increasing function of time. 

It is the purpose of the present paper to consider the evolution of causal bulk viscous dissipative cosmological models 
in the presence of variable gravitational and cosmological constants. By assuming that A is a quadratic function of 
the Hubble parameter H and that the total energy of the Universe is conserved, the gravitational field equations can 
be exactly integrated and two classes of solutions are obtained. The first class describes a barotropic viscous fiuid 
with arbitrary 7, while the second solution corresponds to a stiff cosmological fluid. 

The present paper is organised as follows. The basic equations of the model are written down in Section II. In 
Section III we present the first class of exact solutions. The case of the stiff bulk viscous cosmological fluid is considered 
in Section IV. In Section V we discuss and conclude our results. 

II. GEOMETRY, FIELD EQUATIONS AND CONSEQUENCES 

We consider that the geometry of the Universe filled with a bulk viscous cosmological fluid can be described by a 
spatially flat FRW type metric given by 

ds^ = dt^ - c? (t) [dr^ + (dO^ + sin^ Odf/)'^)] . (1) 

The Einstein gravitational field equations with variable G and A are: 
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R,k - ^g,kR = 8^G(i)T,fc + Ait)g^k■ (2) 
In the following we consider a system of units so that c = 1. 

The bulk viscous effects can be generall y d escribed by means of an effective pressure 11, formally included in 
the effective thermodynamic pressure Pef / [|34| . Then in the comoving frame the energy momentum tensor has the 
components Tq = p,Tl = T| = Tg = —Peff- For the hne element (|^) the Einstein field equations give: 

sf^] ^87rG{t)p + A{t), (3) 



3^ = -47rG(t) {3peff +p) + A{t), (4) 

where a dot denotes the derivative with respect to the time t. 

Taking the covariant derivative of the Einstein field equations with respect to the FRW metric (||) we obtain the 
general conservation law in the presence of variable gravitational and cosmological constants: 

where we have also introduced the Hubble parameter H = a/ a. 

Assuming that the total matter content of the Universe is conserved, = Q, the energy density of the matter 

obeys the usual conservation law: 

p + 3H{p,ff+p) = 0. (6) 

In the presence of the bulk viscous stress 11, the effective thermodynamic pressure term becomes Peff = p + U, 
where p is the thermodynamic pressure of the cosmological fluid. Then the energy conservation equation (^) can be 
split into two independent equations: 

p + 3H{p + p)^ ~3UH, (7) 

87rG/9 + A = 0. (8) 
The causal evolution equation for the bulk viscous pressure is given by 

Tn + n = -3^H-^Tu(^3H + ^-^-^y (9) 

where T is the temperature, ^ the bulk viscosity coefficient and t the relaxation time. 
The growth of the total comoving entropy S(t) over a proper time interval {to,t) is 

I](t)-E(to)--^ / -^dt, (10) 

where /cb is the Boltzmann constant. 

An important observational quantity is the deceleration parameter Q — {jj) ^ ^- The sign of the deceleration 
parameter indicates whether the model inflates or not. The positive sign of q corresponds to "standard" decelerating 
models whereas the negative sign indicates inflation. The deceleration parameter can be expressed as a function of 
the thermodynamic, gravitational and cosmological quantities in the form 

^ AwGit) (p + 3p + 3n)-Aft) 

87rG{t)p + A{t) ■ ^ ' 

Curvature is described by the tensor field Rlj^- It is well-known that if one uses singular behavior of the components 
of the tensor or its derivatives as a criterion for singularites, one gets into trouble since the singular behavior of 
components could be due to singular behavior of the coordinates or tetrad basis rather than the curvature itself. To 
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avoid this problem, one should examine the scalars formed out of curvature. The invariants RijR^^ and RijkiR^^'^'' 
(the Kretschmann scalar) are very useful for the study of the singular behavior of the flat FRW metric: 



R^JR'^ = 12 (iH^ + 3HH^ + H^^ (12) 



R.jkiR'^''^ = 12 (2H^ + 2HH^ + H^j (13) 

In order to close the system of equations (^ and (0)-(^) we have first to give the equation of state for p and specify 
T, T and ^. We shall assume the following laws, p^ : 

p={j-l)p,^ = ap',T = Top^,T=^=ap'-\ (14) 

P 

where 1 < 7 < 2, a > 0, Tq > and s > are constants. Equations (|l^) are standard in the study of bulk viscous 
cosmological models, whereas the equation r — ^/p is a simple procedure to ensure that the speed of viscous pulses 
does not exceed the speed of light. 

With these choices the general solution of the gravitational field equations with variable gravitational and cosmo- 
logical constants still depends on the functional form of G and A. In the present paper we shall fix the mathematical 
form of the cosmological constant assuming that it is a function of the Hubble parameter only and its time dependence 
is: 

A = 3/3^2. (15) 

This ansatz, initially proposed in [ p^ on dimensional ground has been widely used to study decaying vacuum 
cosmological models [Fsl- [^T). 

From Eqs. (^), (||) and ( p5| ) it follows that in this model the gravitational constant G and the energy density p are 
given by: 

G = bH-"^,p = poH", (16) 

where po — ^ and 6 > are constants and n — 2/ {1 — (3) > 0. For physically realistic cosmological models the 
condition < < 1 must hold, in order to assure a time decreasing energy density of the Universe. One must also 
assume that 77 is a decreasing function of the cosmological time. 

With the use of the barotropic equation of state p = (7 — 1) p and Eq. (^, we obtain 

ij + — + 47rGn = 0. (17) 
n 

+ — + 47r6ff-"^n = 0. (18) 



In view of Eq. (|16D, Eq. (|17|) becomes 



With the use of Eqs. (ujm and (18), the causal evolution equation for the bulk viscosity (O) leads to the following 



equation for the Hubble function H: 



H 



J \27/ n \ 2 J n 



where we have denoted ao = (po)^ * /o/-- Eq- (|19|) generalizes the standard evolution equation for causal bulk viscous 
models to the case of variable cosmological and gravitational constants and reduces to it in the case G = const., 
A ~ const.. 

By introducing the transformations 



Eq. ([1^ reduces to 



where we have denoted 



H = y^,r^= / y^dt, (20) 



_ 27[«(l-^)-l] 
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III. THE FIRST CLASS OF SOLUTIONS 



The general solution of the Eq. ( pT| ) depends on the values of the parameter m, combining the effects of the 
thermodynamic equation of state of matter, of bulk viscosity and of the variation of the cosmological constant. Since 
a general solution for arbitrary m is difficult to obtain, we shall limit our study to some particular values of m, for 
which the general solution Eq. can be expressed in an exact analytical form. 

As a first class of solutions we consider the solutions generated by the choice m = 0, corresponding tos = (1 + /3) /2, 
with < (3 < 1 and 1/2 < s < 1. In this case the evolution of the bulk viscous pressure is directly coupled to the 
time variation of the gravitational constant, which is described by the parameter (3. 

Then the general solution of Eq. (^) can be obtained immediately in the form: 



y = c+e"+'" + c^e"-", 



(22) 



where c+ and c_ are constant of integration and we have denoted m± — "° — V°o+^^ 2^ 

Therefore we can express the general solution of the gravitational field equations for a FRW isotropic space-time 
filled with a bulk viscous cosmological fluid in the framework of the full Israel-Stewart-Hiscock causal theory with 
variable cosmological and gravitational constants in the following exact parametric form: 



t-h 



-27/n 



df], 



(23) 



A = 3/3 (c+e'"+'' + c^e'"-")^^/" ,G^b (3/3)"'^/' A-"'^/^ p = piA"/^ 



p = pi (7 - 1) A"/2, e = a {p,r A«"/2^ T = To (pi)^ A^ ,t = a (pi)^"^ At(«-i) , 



n = -ii^-^ (€+6™+" + c-e™-")^'^ 



(24) 
(25) 

(26) 
(27) 

(28) 



37-1 



(29) 



So 



3(l~/3)7ag 

kB^TTbToipo)^ 



(30) 



R.^kiR"^'"^ = 24 (0+6™+" + c.e™-")^'^/" 



n n 



(31) 



where to^ ao and Eq are constants of integration. We have denoted F (77) 



c_m_ e 



c+e + '+c_e 



and PI = ^ 
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IV. THE SECOND CLASS OF SOLUTIONS 



Equation ( pT] ) has also other solutions describing the dynamics of the causal bulk viscous FRW universe for appro- 

dy 
drj 



priate choice of the parameters m and 7. By means of further substitutions v = 1/u and u — ^ the second order 



differential equation (01^) can be transformed into an Abel type differential equation 
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dv 
dy 



y + —y 



+ (3 + ao2/™) «2 



By introducing two new functions A (y) = g(^^/2ti)%^)+3"J'oy../2 and B (y) 
rewrite Eq. in the general form 



9(l-7/2)y 3aoy" 
27 2 



(32) 



allows to 



dv 



dy B {y) 



d Ajy) 
dy B {y) 



By introducing a new variable 



_1 Ajy) 
" V B{yy 



Eq. (p3) can be rewritten in the general form 



^ =aB {y) + A{y)., 



or equivalently 



3 + aoy""/ (m + 1) 



d^ " [9(l-7/2)y/ (27) - 3aoy™+V2] ^ 9 (7/2 - 1) / (27) + 3aoy"V2 ' 



(33) 



(34) 



(35) 



(36) 



Since the energy density of the early universe is supposed to be extremely high, we can consider the extreme limit 
of a stiff cosmological fluid with equilibrium pressure equal to the energy density and then the parameter 7 = 2. In 
the case m = —2, that is 3s = 1 + /3, < /3 < 1 and 1/3 < s < 2/3, Eq. ( |3^ ) becomes a Riccati type differential 
equation: 



dy 2 2 2 2 
da Qfo oao 3 



A particular solution of Eq. ( ^7| ) is given by 

y = -a(aVao + 3/2)"\ 
In the following we denote A (cr) — (cr^/ao + 3/2) ^ . Therefore the general solution of Eq. (^7[) is 



(37) 



(38) 



y = -aA (a) + 



A2 (a) 



Ci-^/A2 (a)e-^da 



(39) 



where Ci is an arbitrary constant of integration. 

The general solution of the gravitational field equations for a flat FRW space-time filled with a bulk viscous 
cosmological fluid and variable G and A can be obtained in the following exact parametric form, with a taken as 
parameter: 



(40) 



H = y2(i-'3) (ct) , a = ao exp 



3an 



y (cr) da 



(41) 



A = 3/3y4(i-« (a) , G = hy-'^^ {a),p^ poy^ {a) , 



P = Poy' (^) , C = a {poY y'' (^) , T = Jo {pof^ y\a),T = a (po)'"' y'^^^'^ ('t) , 



(42) 
(43) 
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9 = 2 (1 - /?) 



(7 tto 



y{a) y^{a) 



1, n = (a) [3y2 (a) - 2ay {a) - 2ao] , 



(44) 



n 


2 


32/2 (<T)-2<Ty((T)-2ao 


p 


~ 3 





(45) 



(l-/3)a3 



2'KhkBOLoTo (po) 



V] a f 2 f 
\\u2 / 2/ (^) [3?/" C'^) - 2<7y (a) - 2ao] exp y (a) da 



da, 



(46) 



RijkiR''"' = 24J/4-8/3 y4 ^ 2 (/3 - 1) y2 (3^2 _ _ ^q) + 2 (/3 - 1)^ (3^^ - ya - ao 



(47) 



where fo, ao and So are constants of integration. 



V. DISCUSSIONS AND FINAL REMARKS 



In the present paper we have considered the evolution of a full causal bulk viscous Universe in the presence of 

variable gravitational and cosmological "constants". Wc have adopted a particular model for the time variation of 
the cosmological constant, assuming that it is a quadratic function of the Hubble parameter H. Consequently, the 
gravitational constant has also a power dependence on the Hubble parameter. In both classes of solutions the evolution 
of bulk viscosity coefficient is related to that of G and A. While the first solution describes a cosmological fluid with 
barotropic equation of state for arbitrary 7, the range of application of the second solution is restricted to a very high 
energy density fluid, obeying the stiff equation of state. 

The causal bulk viscous Universe described by the first class of solutions starts its evolution from a singular state 
with zero value of the scale factor a, a (tg) = 0, and infinite value of the energy density, p{to) — > 00. 

The existence of a singular initial state is also showed by the behavior of the invariant RijkiR^^'^'' , represented in 
Fig.l, which is singular function of time. 
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FIG. 1. Time variation of the invariant RijkiR'^''^ for the first class of solutions for a radiation filled bulk viscous Universe 
(7 = 4/3), for different values of the parameter (3: (5 = 0.1 (solid curve), (5 = 0.3 (dotted curve), P = 0.5 (dashed curve), 
(3 = 0.6 (long dashed curve). We have chosen the integration constants so that c+ = c_ = 1 and assumed that ao = 1. 

The evolution of the Universe is expansionary, with the scale factor an increasing function of time. The Hubble 
parameter, presented in Fig. 2, is a monotonically decreasing function of time, similar to the energy density of the 
cosmological fluid, represented in Fig. 3. 
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FIG. 2. Evolution of the Hubble paxameter H for the first class of solutions for a radiation filled bulk viscous Universe 
(7 = 4/3), for different values of the parameter (3: (3 = 0.1 (solid curve), (3 = 0.3 (dotted curve), (3 = 0.5 (dashed curve), 
(3 = 0.6 (long dashed curve). We have chosen the integration constants so that c+ = c_ = 1 and assumed that ao = 1. 
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FIG. 3. Dynamics of the energy density p for the first class of solutions for a radiation filled bulk viscous Universe (7 = 4/3), 
for different values of the parameter /3: /3 = 0.1 (solid curve), (3 = 0.3 (dotted curve), /3 = 0.5 (dashed curve), (3 = 0.6 (long 
dashed curve). We have chosen the integration constants so that c+ = c- = 1 and assumed that ao = 1. 

The bulk viscous pressure 11, shown in Fig. 4, satisfies the condition 11 < only for time intervals greater than an 
initial value ti. 




FIG. 4. Variation of the bulk viscous pressure 11 of the cosmological fluid for the first class of solutions for a radiation filled 
bulk viscous Universe (7 = 4/3), for different values of the parameter (3: (3 = 0.1 (solid curve), /3 = 0.3 (dotted curve), /3 = 0.5 
(dashed curve), (3 = 0.6 (long dashed curve). We have chosen the integration constants so that c+ = c- = 1 and assumed that 
ao = 1. 

Hence this model can describe the dynamics of the causal bulk viscous Universe with variable gravitational and 

cosmological constants only for a finite time interval. 

The deceleration parameter q, represented in Fig. 5, shows an initial non-inflationary evolution, with q > 0, but in 
the large time limit the Universe ends in an inflationary epoch. 
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FIG. 5. Dynamics of the deceleration parameter q for the first class of solutions for a radiation filled bulk viscous Universe 
(7 = 4/3), for different values of the parameter (3: (3 = 0.1 (solid curve), /3 = 0.3 (dotted curve), (3 = 0.5 (dashed curve), 
/? = 0.6 (long dashed curve). We have chosen the integration constants so that c+ = c- = 1 and assumed that ao = 1. 

During the viscous effects dominated era a large amount of comoving entropy is produced. Fig. 6 shows the time 
variation of the entropy. 
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FIG. 6. Time variation of the comoving entropy E for the first class of solutions for a radiation filled bulk viscous Universe 

(7 — 4/3), for different values of the parameter 13: j3 = 0.1 (solid curve), /? = 0.3 (dotted curve), f3 — 0.5 (dashed curve), 
(3 — 0.6 (long dashed curve). We have chosen the integration constants so that c+ = C- = 1 and assumed that ao = 1. 

The time variation of the cosmological and gravitational constants is represented in Figs. 7 and 8. 
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FIG. 7. Time variation of the cosmological constant A{t) for the first class of solutions for a radiation filled bulk viscous 
Universe (7 = 4/3), for different values of the parameter (3: (3 = 0.1 (solid curve), (3 = 0.3 (dotted curve), (3 = 0.5 (dashed 
curve), (3 = 0.6 (long dashed curve). We have chosen the integration constants so that c+ = c- = 1 and assumed that ao = 1- 
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FIG. 8. Time variation of the gravitational constant G{t) for the first class of solutions for a radiation filled bulk viscous 
Universe (7 = 4/3), for different values of the parameter (3: (3 = 0.1 (solid curve), /3 = 0.3 (dotted curve), /3 = 0.5 (dashed 
curve), /3 = 0.6 (long dashed curve). We have chosen the integration constants so that c+ = c- = 1 and assumed that ao = 1. 

The cosmological constant is a decreasing function of time, while the gravitational constant G tends in the large time 

to a constant value. But in this model the ratio of the bulk viscous pressure and of the thcrmodynamical presssurc 
is greater than one during the inflationary period. Consequently during this period the model is not consistent 
thermodynamically. 

The evolution of the stiff cosmological fluid filled bulk viscous Universe, described by the second class of solutions, 
also starts from a singular state, with a {to) = 0. The invariant RijkiR^-'^^ , represented in Fig. 9, is also singular in 
the large time limit. 
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FIG. 9. Time variation of the invariant RijkiFC^'^'' for the second class class of solutions describing a stiff fluid flUcd bulk 
viscous Universe (7 = 2), for different values of the parameter /3: /J = 0.1 (solid curve), (3 = 0.3 (dotted curve), (3 = 0.5 (dashed 
curve), (3 = 0.6 (long dashed curve). We assumed that ao = 1 and we have chosen the arbitrary integration constant Ci = 1/2. 

The evolution of the Universe is expansionary, with the scale factor, represented in Fig. 10, an increasing function 
of time. 
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FIG. 10. Time evolution of the scale factor a for the second class class of solutions describing a stiff fluid filled bulk viscous 
Universe (7 = 2), for different values of the parameter /3: /3 = 0.1 (solid curve), /3 = 0.3 (dotted curve), /3 = 0.5 (dashed curve), 
(3 = 0.6 (long dashed curve). We assumed that ao = 1 and we have chosen the arbitrary integration constant Ci = 1/2. 

The Hubble parameter and the energy density, represented in Fig. 11, are monotonically decreasing functions of 
time. 
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FIG. 11. Dynamics of the energy density p for the second class class of solutions describing a stiff fluid filled bulk viscous 
Universe (7 = 2), for different values of the parameter (3: (3 = 0.1 (solid curve), /3 = 0.3 (dotted curve), /? = 0.5 (dashed curve), 
/? = 0.6 (long dashed curve). We assumed that ao = 1 and we have chosen the arbitrary integration constant Ci = 1/2. 

The bulk viscous pressure 11 satisfies the condition 11 < for all time intervals and for all values of the parameter, 
as can be seen from Fig. 12. 
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FIG. 12. Time variation of the bulk viscous pressure II for the second class class of solutions describing a stiff fluid filled 
bulk viscous Universe (7 = 2), for different values of the parameter (3: f3 = 0.1 (solid curve), f3 = 0.3 (dotted curve), /3 = 0.5 
(dashed curve), (3 = 0.6 (long dashed curve). We assumed that ao = 1 and we have chosen the arbitrary integration constant 
Ci = 1/2. 

The deceleration parameter q, represented in Fig. 13, satisfies in this model the condition g > for all times and 
in the acceptable range of the physical parameters. 
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FIG. 13. Time dependence of the deceleration parameter q in the second class class of solutions describing a stiff fluid filled 
bulk viscous Universe (7 = 2), for different values of the parameter /3: f3 = 0.1 (solid curve), /? = 0.3 (dotted curve), /? = 0.5 
(dashed curve), /3 = 0.6 (long dashed curve). We assumed that ao = 1 and we have chosen the arbitrary integration constant 
Ci = 1/2. 

Therefore the evolution of the high density causal bulk viscous Universe is non-inflationary. Consequently, the 
model is consistent therniodynamically satisfying the condition of the smallness of the bulk viscous pressure. The 
time evolutions of the cosmological constant represented in Figs. 14 is similar to the behavior of this quantity in the 
first model, with the cosmological constant a decreasing function of time, while in the large time limit the gravitational 
constant, showed in Fig. 15, is a slowly increasing function of time. 
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0.05 0.1 0.15 0.2 0.25 0.3 
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FIG. 14. Time variatiou of the cosmological constant A for the second class class of solutions describing a stiff fluid filled 
bulk viscous Universe (7 = 2), for different values of the parameter (3: (3 = 0.1 (solid curve), (3 = 0.3 (dotted curve), (3 = 0.5 
(dashed curve), /J = 0.6 (long dashed curve). We assumed that ao = 1 and we have chosen the arbitrary integration constant 
Ci = 1/2. 
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FIG. 15. Evolution of the gravitational constant G(t) for the second class class of solutions describing a stiff fluid filled bulk 
viscous Universe (7 = 2), for different values of the parameter /3: (3 = 0.1 (solid curve), (3 = 0.3 (dotted curve), (3 = 0.5 (dashed 
curve), (3 = 0.6 (long dashed curve). We assumed that ao = 1 and we have chosen the arbitrary integration constant Ci = 1/2. 
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The possibility that the cosmological constant and the gravitational coupling are not real constants is an intriguing 
possibility, which has intensively been investigated in the physical literature. It is a very plausible hypothesis that 
these effects were much stronger in the early Universe, when dissipative effects also played an important role in the 
dynamics of the cosmological fluid. Hence the solutions obtained in the present paper could give an appropriate 
description of the early period of our Universe. 
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